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CROSSFLOW IN TWO-DIMENSIONAL ASYMMETRIC NOZZLES 


Daniel I. Sebacher and Louise P. Lee 
Langley Research Center 

SUMMARY 

An experimental investigation of the crossflow effects in three contoured, two- 
dimensional asymmetric nozzles is described and results are compared with theoretical 
predictions of nozzle flow by use of an inviscid method of characteristics solution and 
two-dimensional turbulent boundary- layer calculations. The nozzles varied in maximum 
expansion angle from 15° 17’ to 29°20' but all were designed for the same area ratio to 
give an exit Mach number of 3.49 by use of expanding air at a stagnation temperature of 
290 K. The results of this study show that crossflow takes place in the boundary layer of 
all the test nozzles. Exit pitot -pressure surveys and sideplate oil drop flows illustrate 
the trade-off between increased crossflow driving force in the larger turning angle nozzle 
to the thicker boundary -layer buildup in the longer, smaller turning angle nozzle. Static 
wall -pressure measurements were found to agree well with the inviscid method of char- 
acteristics theory when it is assumed that no crossflow and a relatively undisturbed flow 
core suitable for wind-tunnel testing exists for the three nozzles. 

INTRODUCTION 


Test facilities (ref. 1) tailored to investigate rectangular scramjet engine configura- 
tions for the integrated scramjet engine -airframe concept (ref. 2) will probably utilize a 
rectangular nozzle exit which is generally obtained by some type of two-dimensional 
expansion. Two-dimensional asymmetric nozzle configurations consisting of a contoured 
wall, a flat bottom plate, and flat sidewalls have been considered for use in both the inte- 
grated scramjet engine -airframe concept and the wind-tunnel test facility and are the 
subject of this report. One problem associated with this type of nozzle is the development 
of lateral pressure gradients on the nozzle sidewalls induced by streamline curvature 
which can result in boundary -layer crossflow along the sidewalls in a direction normal to 
the potential flow streamlines. These lateral pressure gradient effects will appear as 
lower pressure along the contoured wall compared with the pressure at the same station 
along the bottom plate in the region of the initial circular -arc expansion and as higher 
pressures along the contoured wall when the flow is turned back parallel downstream of 
the maximum turning angle. The nozzles designed for this study will also contain a region 



of parallel flow where design pressure has been achieved first along the bottom plate and 
then through the flow core to the nozzle exit on the contoured wall. In this region of uni- 
form pressure, the previous driving force for crossflow vanishes, but the difference in 
streamline direction between the flow core and the boundary layer will further modify the 
boundary -layer flow by viscous interaction. 

Since the use of the two-dimensional asymmetric nozzles has been limited and 
generally confined to small expansion angles, the effects of crossflow in these nozzles have 
largely been ignored. A related investigation of secondary flows in the convergent region 
of planar nozzles is reported in reference 3, and an analysis of crossflow in laminar flow 
with boundary-layer control using suction is given in reference 4. A novel study of the use 
of fences to decrease crossflow effects on a nozzle sidewall is found in reference 5. 

To determine the effects of crossflow as a function of maximum nozzle turning angle, 
the flow characteristics of three contoured two-dimensional asymmetric nozzles are inves- 
tigated experimentally and the results are compared with nozzle flow predictions by using 
an inviscid method of characteristics solution and two-dimensional boundary -layer theory. 
Any significant distortion of the flow core due to crossflow is of particular interest because 
of the possible application of nozzles of this type to integrated engine-airframe, hypersonic 
aircraft configurations. (See ref. 2.) 
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SYMBOLS 

throat height (see fig. 1(b)) 

Reynolds number (based on flow speed at nozzle throat and throat height) 

static pressure on nozzle bottom plate 

static pressure on nozzle contoured wall 

stagnation chamber pressure 

impact pressure 

static pressure on sidewalls 

nozzle coordinates (see fig. 1(a)) 

nozzle maximum expansion angle 
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Abbreviations: 


2-D two-dimensional 

MOC method of characteristics 

APPARATUS AND METHOD 
Nozzles 

A preliminary design analysis, using the method of characteristics solution, indi- 
cated that for the flow conditions available for this experiment, a suitable range of the 
sidewall pressure gradients causing crossflow in two-dimensional asymmetric nozzles 
could be achieved with # max between 15° and 29°. Therefore, three contoured nozzles 
with nominal values of 6 of 15°17', 23°23', and 29°20' were selected for investiga- 
tion as shown in figure 1. The three nozzles were identical upstream of the throats and 
all were designed to expand 290 K air to a free-stream Mach number of 3.49. The length 
of each nozzle wall was terminated at the location that theoretical calculations indicated 
parallel flow would be achieved; this termination resulted in a decreasing nozzle length as 
the maximum expansion angle was increased. The length of the nozzles were 42.10 cm, 
31.10 cm, and 26.62 cm for nozzles with maximum expansion angles of 15°17’, 23°23', 
and 29°20’, respectively. For the three nozzles, the width was 7.62 cm, the throat height 
was 0.847 cm, and the exit height was 6.10 cm. The stagnation chamber housed a perma- 
nent flow straightener, as shown in figure 1(b). The nozzle sidewalls were sealed at the 
contact surfaces of the contoured top plate and the flat bottom plate. 

Instrumentation 

The strain-gage pressure transducers used in this experiment have a rated measure- 
ment accuracy of ±0.5 percent. The output from the transducers is recorded on tape and 
converted to a digitized readout by using a calibration tape; therefore, the usual errors 
involved in reading a visual recording are not applicable. On the other hand, there are 
errors involved in the recording readout system due to instrument drift and systematic 
calibration inaccuracies which have been estimated to be ±2 percent by comparing the 
system readout with a number of direct measurements. The random error of the record- 
ing system due to noise is given as ±0.2 percent by the manufacturer. The maximum error 
for the pressure ratios presented in this report is ±5.4 percent. 

The impact pressure probe used to survey the boundary layer and flow core at the 
nozzle exit plane had a circular cross section with an outside diameter of 0. 1 cm. The 
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probe was 5.23 cm long and was attached to a wedge-shaped strut of the same thickness 
which, in turn, was remotely controlled by a motor -driven screw which permitted travel 
in a direction normal to nozzle walls. The impact -pressure measurements were recorded 
on digital tape and sampled every 0.1 second and thus gave a data point every 0.044 cm in 
the y-direction. A traverse was taken every 0.5 cm in the z-direction from the sidewall 
to 0.5 cm past the center line of the nozzle. 

The static wall pressure orifices have a diameter of 0.15 cm and are connected to 
the pressure transducers by a short piece of stainless-steel tubing of the same inside 
diameter. The coordinates of the static wall -pressure orifices are found in tables I 
and n. The wall static pressures and impact pressures were measured during the same 
runs and recorded on tape along with the stagnation chamber pressure, probe position, and 
airflow rate. 

The oil -flow studies consisted of splattering small quantities of a mixture of lamp- 
black and silicone oil over the area of the sideplates before each run, and observing the 
flow pattern of the mixture after a quick release valve initiated the airflow. A plexiglass 
sideplate was used on one side of the nozzle so that the oil flow could be observed during 
the run. Several seconds of flow time were required to obtain a useful flow pattern which 
could be photographed and once the flow pattern was set, no drift was observed when left 
to stand over a period of days if a proper mixture of lampblack and oil was initially 
applied. When a thin coat of silicone oil was rubbed into the sidewall prior to splattering 
the oil -flow mixture, the quality of the flow patterns was improved. 

Test Conditions and Procedures 

All three nozzles expanded 290 K dried air, with an H 2 O mole fraction less than 
2 x 10“ 4 , to a free- stream Mach number of 3.49. The nozzle flow was exhausted into an 
evacuated chamber maintained at 5 mm Hg. Stagnation pressure varied from 220 to 
372 kN/m^ and resulted in a Reynolds number range of 2.75 x 10^ to 4.64 x 10^ based on 
flow conditions at the nozzle throat and throat height. These conditions along with the 
physical size of the nozzles were largely limited by the maximum mass flow of dry air 
available which was 0.68 kg/sec. 

METHOD OF CALCULATION 
Nozzle Contouring 

The numerical technique for determining the contours for the two-dimensional 
asymmetric nozzles used the method of characteristics for a rotational gas mixture. The 
stagnation enthalpy is allowed to vary normal to the streamlines. (See refs. 6 and 7.) 

Since only expanding 290 K air was considered, this analysis assumes each component of 
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the gas to be thermally perfect. Properties for the chemical species of air as a function 
of temperature, as needed for this method, are presented in reference 8. Initial uniform 
profiles of velocity, temperature, pressure, flow inclination, and mass fraction at the 
nozzle throat for starting the contouring program are also calculated by assuming a ther- 
mally perfect gas. (See ref. 9.) 

When using the method of characteristics solution for computing the three nozzle 
contours shown in figure 1(b), the flow was initially expanded through a circular -arc 
contour of specified radius to a point of maximum turning angle from which the down- 
running characteristic line yields the design pressure on the nozzle axis. For the flow 
conditions and exit Mach number of these nozzles, the maximum turning angle is a strong 
function of the initial circular arc up to a 0 max of about 29°, after which the maximum 
turning angle becomes insensitive to decreasing the initial circular arc. 

The method of characteristics solution will also predict the inviscid static -pres sure 
distribution on the nozzle walls and throughout the flow field. The results of these calcu- 
lations were not only used as a comparison with the measured static pressures, but were 
also needed as input to the boundary -layer calculations. The analytical program used to 
compute the nozzle contours had a number of other computational options available but 
only those needed in this study have been described. Proper grid spacing was determined 
to insure accuracy in these solutions and the program’s iteration procedure usually con- 
verged in three iterations. Coordinates for the three nozzles investigated are found in 
table HI. 


Boundary -Layer Calculations 

Solutions of two-dimensional laminar and turbulent boundary -layer equations 
for multicomponent nonreacting gases were used in the boundary -layer analysis. The 
boundary -layer equations are transformed by the Levy- Lees transformation equations. 
These solutions (ref. 10) used an implicit finite-difference scheme for two-dimensional 
nozzles. The energy equation is written in terms of the stagnation enthalpy and the 
Reynolds' shear, and conductivity terms are given by an eddy viscosity and the turbu- 
lent Prandtl number. 

The turbulent analysis employs the eddy viscosity law derived by Reichardt in the 
region of the boundary layer near the wall and follows Clauser's work using Klebanoff’s 
inter mittency factor in the outer part of the boundary layer. Reference 11 presents the 
system of equations describing the fluid flow in a boundary layer and a FORTRAN pro- 
gram used to solve these equations. This program has the capability of handling two- 
dimensional laminar or turbulent boundary -layer flows and the transition from laminar to 
turbulent flow. The computational step size along the wall is varied throughout the nozzle 
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depending on the pressure gradient at each point. The boundary -layer wall pressure 
distribution was determined by using the inviscid method of characteristics solution. 

RESULTS AND DISCUSSION 

Nozzle Differential Pressure 

The nondimensional differential pressures between the contoured wall (P cw ) and the 
bottom plate ^p^ of the asymmetric nozzles computed by the inviscid theory are shown 
in figure 2 as a function of axial distance. These pressure differentials for the three 
nozzles illustrate the crossflow driving force in terms of pressure gradient variations 
along the nozzle sidewalls. The initial dip in the pressure differential takes place because 
the flow first expands more rapidly along the circular -arc contour than along the bottom 
plate. Once the wall curvature reverses and begins to turn the flow back parallel to the 
nozzle bottom plate, the pressure on the contoured wall becomes greater than the pressure 
on the bottom plate. The reversal of the streamline curvature therefore causes the 
pressure differential to increase to zero, then go positive, and again decrease to zero as 
the nozzle exit is approached. Both the position of the pressure differential peak and its 
magnitude should affect the degree of crossflow in these nozzles. The maximum pressure 
differential not only takes place at a point of greater pressure in the higher turning angle 
nozzles but also farther upstream in the critical region just downstream of the throat. 

This maximum pressure differential decreases in magnitude and its locus moves down- 
stream as the turning angle is decreased. The nozzles designed for this study also con- 
tain a region of parallel flow where design pressure has been achieved first along the 
bottom plate and through the flow core to the nozzle exit on the contoured wall as indicated 
by the area downstream of the parallel flow characteristic lines in figure 1(b). In this 
region of uniform pressure the previous driving force for crossflow vanishes. 

Oil Flows 

Experimental verification of the boundary -layer crossflow normal to the inviscid 
flow streamlines is shown in the oil -flow studies presented in figure 3. Any upward 
crossflow from the bottom plate toward the contoured wall due to the initial negative pres- 
sure gradient indicated in figure 2 is not readily discernible for any of the nozzles in 
figure 3. On the other hand, the downward crossflow from the contoured wall to the bottom 
plate is obvious in these photographs for all three nozzles because of the large broad pos- 
itive gradients shown in figure 2. The streamlines start to become parallel to the bottom 
plate as the exit is approached. In this region of uniform pressure the crossflow driving 
force vanishes, but the difference in streamline direction between the flow core and the 
boundary layer will further modify the boundary -layer flow direction by viscous inter- 
action. Additional crossflow effects have been observed in the test sections of the 
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two-dimensional nozzles with constant cross-section extensions (refs. 12 and 13) but the 
flow region past the nozzle exit is not considered in this study. A typical predicted 
streamline with no crossflow has been included in figure 3 to indicate the direction of the 
core flow. 


Wall -Pressure Measurements 

Wall -pressure measurements were taken along the center line of the contoured plate 
and the bottom plate, and also at various locations on a grid inscribed on the sideplate for 
each nozzle configuration. (See tables I and n.) These data are presented as the pres- 
sure difference between the contoured and bottom plates as a function of axial distance in 
figure 4 and as sidewall isobars in figure 5. The inviscid two-dimensional pressure dif- 
ferences and isobars as predicted by the method of characteristics solutions are also 
shown in these figures. The measured pressure differences in figure 4 generally agree 
with the inviscid theory which does not account for boundary -layer crossflow effects even 
though crossflow is known to occur along the sidewalls. The pressure differential meas- 
ured on the center line of the top and bottom plates is generally less than the pressure 
differential measured on the top and bottom of the sidewalls. 

Sidewall pressure measurements similarly agreed with the static pressures pre- 
dicted by inviscid theory as seen in figure 5. In figure 5 the theory and data are presented 
in isobar form with the mirror image of the asymmetric nozzle comparing the experimen- 
tal data with the theoretical predictions. The contours drawn through the experimental 
data by using the theory as a guide were drawn as smooth curves because of the uncertainty 
of determining the pressure gradients between the pressure orifices. 

The agreement between the inviscid predictions and the measured wall pressures 
indicate that even when crossflow occurs in the boundary layer, the inviscid pressure gra- 
dients driving the crossflow could be used in a crossflow theoretical model, at least for 
the conditions of this experiment. The discrepancy that does exist between the predicted 
and measured wall-pressure distributions results from no allowance for the boundary- 
layer growth on the flat walls in the design of the nozzles. A boundary -layer correction 
was applied to the contoured wall. 

Impact- Pressure Surveys 

Impact-pressure traverses were made in the y-direction at constant values of the 
z -coordinate at the exit plane of each nozzle. A large number of traverses were measured 
of which figures 6 and 7 are typical. The numerous points measured in the exit plane of 
each nozzle were then plotted as a function of y and z and isobars were contoured 
through the data as illustrated in figure 8. These impact pressure traverses and isobar 
contours suggest the effects of crossflow in the boundary layer by the existence of 
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total pressure losses along the sidewalls. In the traverse of figure 6, the two-dimensional 
turbulent boundary -layer analysis predicts the total -pressure distribution very well with 
the exception of the outer region of the boundary layer on the flat bottom plate. The tra- 
verse shown in figure 7 indicates the pressure loss due to crossflow very near the sidewall. 
The isobars of figure 8 show lines of constant values of impact pressure ratio for half of 
the exit region of each nozzle. The indication of secondary flow down the sidewall is 
clearly discernible. 

For the asymmetric nozzles studied in this experiment, the crossflow effects appear 
to be generally contained in the boundary layer. An interpretation of the crossflow effects 
in the corner regions of the nozzle will not be attempted in this analysis. The contours of 
figure 8 also show that the crossflow effects at the exit are no more severe for the higher 
turning angle nozzle even though the driving force is larger. This condition is due to the 
compensating effect between increased crossflow driving force in the large turning angle 
nozzle as compared with the thicker boundary -layer buildup in the longer, but smaller 
turning angle nozzle. Since the crossflow effects are largely contained in the boundary 
layer, a reasonable undisturbed flow core was measured for the three test nozzles. A 
comparison of the flow core found in these asymmetric nozzles with the flow core in a full 
two-dimensional nozzle may be made by using the data of references 5 and 12. In the sym- 
metrical two-dimensional nozzle, the test flow core is reduced by the distorted boundary- 
layer growth along the center line of the sidewalls due to boundary -layer crossflow from 
both contoured walls toward the sidewall center line. 

Reynolds Number Effects 

All the data shown up to figure 8 were taken at a Reynolds number of 4.64 x 10®, 
based on the flow speed at the nozzle throat and the throat height. The experiment was 
then repeated for Reynolds numbers of 3.70 x 10 5 and 2.75 x 10® by reducing the nozzle 
stagnation pressure from 372 to 220 kN/m 2 and maintaining an exit pressure to allow full 
expansion in the nozzle. Figures 9 and 10 illustrate the typical effects of Reynolds num- 
ber on the total-pressure traverse at the exit of the 23°23' nozzle. The traverse at 
z = 0.107 mm from the sidewall in figure 9 for Reynolds numbers of 4.64 x 10® and 
2.75 x 10® indicate a slightly greater total-pressure loss in the bottom plate region for 
the lower Reynolds number. The traverse at z = 0.5 mm from the sidewall in figure 10 
which is inside the sidewall boundary layer indicates an even greater total -pres sure loss 
at the lower Reynolds number. The increase in pressure loss is expected since the wall 
boundary -layer thickness increases with a decreasing Reynolds number. The effect of 
crossflow on this pressure loss is difficult to assess. Over the limited range of Reynolds 
number covered in this experiment, the effects are very small and the rest of the impact 
data is similar to that shown in figures 9 and 10. The static wall-pressure data show no 
Reynolds number effects as expected. 
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SUMMARY OF RESULTS 


The crossflow effects in three contoured, two-dimensional, asymmetric nozzles have 
been experimentally investigated and the results are compared with theoretical predictions 
of nozzle flow by using an inviscid method of characteristics solution and two-dimensional 
turbulent boundary-layer calculations. The results of this study show that 

1. Crossflow takes place in the three test nozzles including a nozzle with a maximum 
turning angle of 15°17'. The existence of crossflow on the sidewalls was verified by oil- 
flow studies. 

2. Crossflow effects are mainly confined to the boundary -layer region of the two- 
dimensional aysmmetric nozzles. The boundary -layer distortions are evaluated by 
comparing the measured impact pressure variations with the predicted two-dimensional 
turbulent boundary -layer profiles. For a nozzle design, where crossflow is undesirable, 
there is a trade-off between increased crossflow driving force in the larger turning angle 
nozzle and the thicker boundary -layer buildup in the longer, smaller turning angle nozzle. 

3. Static wall pressure measurements agree well with the predictions of the inviscid 
two-dimensional theory even though crossflow occurs and is not accounted for in the 
calculations. 

4. A relatively undisturbed flow core suitable for wind-tunnel test facilities and 
rectangular engine nozzle configurations was found for the three test nozzles. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, Va. 23665 
May 30, 1975 
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TABLE I.- SIDEWALL PRESSURE ORIFICE COORDINATES 
















TABLE II.- AXIAL DISTANCES FOR PRESSURE ORIFICES 
ALONG CENTER LINE OF CONTOURED PLATE 
AND FLAT BOTTOM PLATE 


— 
x, cm 

x, cm 

x, cm 

= 15°17' nozzle 

IIlctA 

0 = 23°23' nozzle 

IIldA 

0 max = 29 ° 20 ' nozzle 

3.38 

1.73 

0.38 

7.29 

3.36 

1.30 

10.70 

5.02 

2.60 

14.10 

6.78 

3.85 

17.50 

10.10 

5.10 

24.28 

13.48 

6.80 

31.08 

16.88 

10.18 

37.80 

23.63 

13.54 


30.45 

17.16 



24.03 












TABLE m.- NOZZLE COORDINATES 


(a) Approach (b) Expansion 


x, cm 

y, cm 

x, cm 

y, cm 

x, cm 

y, cm 

0 max = 1 5° 17 1 nozzle 

0 max = 23°23’ nozzle 


0.2540 

0.85598 

0.2540 

0.85598 

0.2540 

0.89916 

.5080 

.86360 

.5080 

.86614 

.5080 

1.02616 

.7620 

.86614 

.7620 

.88646 

.7620 

1.17094 

1.0160 

.87122 

1.0160 

.91694 

1.0160 

1.31572 

1.2700 

.87884 

1.2700 

.95250 

1.2700 

1.46050 

1.5240 

.88900 

1.5240 

.99822 

1.5240 

1.59512 

1.7780 

.90170 

1.7780 

1.05410 

1.7780 

1.73228 

2.0320 

.91440 

2.0320 

1.11506 

2.0320 

1.85928 

2.2860 

.92710 

2.2860 

1.18618 

2.2860 

1.99390 

2.5400 

.94234 

2.5400 

1.25984 

2.5400 

2.11582 * 

3.0480 

.97536 

3.0480 

1.44526 

3.0480 

2.34950 

3.5560 

1.01854 

3.5560 

1.65354 

3.5560 

2.57302 

4.0640 

1.06680 

4.0640 

1.86690 

4.0640 

2.78130 

4.5720 

1.12014 

4.5720 

2.07264 

4.5720 

2.98450 

5.0800 

1.17348 

5.0800 

2.27330 

5.0800 

3.16992 

6.3500 

1.34620 

6.3500 

2.72288 

5.5880 

3.33756 

7.6200 

1.57226 

7.6200 

3.13436 

6.0960 

3.50012 

8.8900 

1.84404 

8.8900 

3.49504 

6.6040 

3.65252 

10.1600 

2.15646 

10.1600 

3.82778 

7.1120 

3.79984 

11.4300 

2.48920 

11.4300 

4.13004 

7.6200 

3.93700 

12.7000 

2.82956 

12.7000 

4.39420 

8.8900 

4.26974 

13.9700 

3.15468 

13.9700 

4.63042 

10.1600 

4.53898 

15.2400 

3.46202 

15.2400 

4,84378 

11.4300 

4.78282 

16.5100 

3.73126 

16.5100 

5.03936 

12.7000 

5.00634 

17.7800 

3.98780 

17.7800 

5.20700 

13.9700 

5.19938 

19.0500 

4.23672 

19.0500 

5.36194 

15.2400 

5.36956 

20.3200 

4.45516 

20.3200 

5.49910 

16.5100 

5.51180 

21.5900 

4.65328 

21.5900 

5.61594 

17.7800 

5.63372 

22.8600 

4.83870 

22.8600 

5.72262 

19.0500 

5.75056 

24.1300 

5.01650 

24.1300 

5.81406 

20.3200 

5.84200 

25.4000 

5.18160 

25.4000 

5.90042 

21.5900 

5.91820 

26.6700 

5.32892 

26.6700 

5.96646 
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(b) Side view of three asymmetric nozzles investigated. 
Figure 1.- Asymmetric nozzles. 
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Figure 2.- Predicted pressure differences between top contoured plate and asymmetric 
nozzle bottom plate plotted against axial distance by using 2-D inviscid method of 
characteristic solutions. 
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Figure 3.- Oil -flow visualization showing crossflow effects. 
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Figure 3.- Continued. 


Predicted streamline with no crossflow 
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( c ) Lov = 29°29' nozzle, 
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Figure 3.- Concluded. 
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On center line of top and bottom plates 
On top and bottom of sidewalls 
Inviscid 2-D theory (MOC) 
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Figure 4.- Continued. 
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Figure 5. - Measured and predicted pressure contours on side plates. 
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Figure 5.- Continued. 
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Figure 6. - Typical measured and predicted impact pressure ratio profiles at 
exit plane against vertical station at z = 0.107 mm from side plate. 
e ma *= 23 ° 2 3’; =4.64 x 105. 







Figure 7.- Typical measured impact pressure ratio profiles at exit plane plotted 
against vertical station at z *■ 0.5 mm from side plate. 0 m ax = 23°23'; 

Nrc = 4.64 x 10 5 . 
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Figure 10. - Typical measured impact pressure ratio profiles at exit plane plotted 
against vertical station at z = 0.5 mm from side plate. 0 max = 23°23'. 
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